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Adiabaticity and the Fate of Non-Gaussianities: The Trispectrum and Beyond 
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Extending the analysis of [J] beyond the bispectrum, we explore the superhorizon generation of 
local non-gaussianities and their subsequent approach to adiabaticity. Working with a class of two 
field models of inflation with potentials amenable to treatment with the 5N formalism we find that, 
as is the case for f^i^'^^, the local trispectrum parameters tnl and (7nl are exponentially driven 
toward values which are slow roll suppressed if the fiuctuations are driven into an adiabatic mode 
by a phase of effectively single field infiation. We argue that general considerations should ensure 
that a similar behavior will hold for the local forms of higher point correlations as well. 



I. INTRODUCTION 

The observed power spectrum of temperature fluctua- 
tions in the cosmic microwave background radiation 
has provided dramatic evidence for inflation as the origin 
of cosmological fluctuations [^-Q. However, probing the 
inflationary paradigm in detail requires that we study 
more than just the two-point correlation function; we 
must also examine the non-gaussian contributions to the 
primordial fluctuations. 

If we model inflation as being driven by canoni- 
cally normalized slowly rolling scalar fields, these fields 
will generally have gaussian fluctuations at horizon exit 
d, For single field models, local form non-gaussian 
contributions to the curvature perturbations are slow-roll 
suppressed [111, [l^l- With multiple fields this may not 
be the case. The continued evolution of the curvature 
perturbation outside the horizon can mix initially gaus- 
sian scalar fluctuations to produce a large non-gaussian 
component in the spectrum. In particular, it has been 
argued that the quantity /nl**' (characterizing the lo- 
cal form bispectrum) can become large in a number of 
multifield models, through the evolution of the curvature 
perturbation after horizon exit; for a non-exhaustive list 
of examples see [l3l - [27j , or for reviews see [1^ [2^ . 

Since the curvature perturbation continues to evolve 
outside the horizon in multifield infiation, to make pre- 
dictions we should follow inflationary observables until 
they become conserved - a process that occurs once cos- 
mological fluctuations enter a purely adiabatic mode fso'^ 
Issj - or until they are observed. It is important to stress 
the difference between conservation of the curvature per- 
turbation, and its intermittent constancy during some 
periods of the early universe. By conservation we mean 
that for a particular mode of wave number k the cur- 
vature perturbation remains constant while the mode is 
outside the horizon ^ <^ H. The curvature perturba- 
tion may become constant without being conserved, as 
is the case when multiple field inflation takes place along 
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a straight trajectory in held space, or when there exist 
two decoupled fluids with the same equation of state. 
The important distinction in these cases is that the cur- 
vature perturbation may evolve during a later period if 
non-adiabatic fluctuations are present when the inflaton 
traverses a turn, or when one decoupled fluid changes its 
equation of state, say by becoming non-relativistic. We 
will focus in this paper on the conservation of the curva- 
ture perturbation, as this guarantees that correlations are 
insensitive to the subsequent evolution of the universe. 
In a recent paper [l| we explored the effect of the ap- 
proach to adiabaticity on non-gaussianities produced by 
superhorizon evolution in multifield inflation. In the lim- 
ited, but illustrative, class of models where isocurvature 
is eliminated by a short period of single field inflation be- 
fore reheating, we found that the approach to adiabatic- 
ity damped out the leading contribution to /nl'^'- Our 
goal in this paper is to carry the analysis further and con- 
sider the evolution of higher order correlation functions 
in the same class of models. 

In addition to allowing for superhorizon evolution of 
the curvature perturbation, non-adiabatic fluctuations 
which persist through the radiation-dominated era leave 
observable effects on the cosmic microwave background 
[33 - l37| . Current observations are consistent with a 
purely adiabatic power spectrum, though they allow for 
a small contribution from non-adibatic fluctuations 2] . 

Although higher order correlations come with more sig- 
nificant observational challenges, there are weak bounds 
on some limits of the trispectrum, bounds which should 
improve considerably with future experiments. From the 
5 year WMAP data Smidt et al. jsl] have found: 
-7.4 < 5nl/10^ < 8.2 and -0.6 < tnl/IO^ < 3.3, 
and Fergusson, Regan, and Shellard found —5.4 < 
<?nl/10^ < 8.6. (7nl and tnl characterize the tris pec trum 
in a particular "local- form" parametrization (iol. l4l| , the 
details of which we will discuss below. Upcoming exper- 
iments will do significantly better - the^Planck satellite 
should be able to achieve Atnl 560 [41] , and probes of 
large scale structure may be able to provide even tighter 
bounds [i^-lil. 

One might, with some apparent justification, wonder 
why we should bother with the trispectrum. After all 
the bispectrum is guaranteed to vanish in the presence of 
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purely gaussian fluctuations making it a sensitive probe 
of non-gaussianity and, given the challenges in measur- 
ing higher order correlations, one might believe it un- 
likely that we would gain much by the effort. This is 
not, however, the case. Because of the existence of con- 
sistency relations between the bispectrum and trispec- 
trum, measurement of tnl can be a sharper probe of 
non-gaussianity than /nl in some models (if /nl ^ 50) 
and even holds the possibility of ruling out large classes 
of multifield models. For this reason, we would like to 
explore the evolution of the trispectrum in multifield in- 
flationary models as adiabaticity is approached. 

As in our previous work we use the SN formalism to 
follow the evolution of perturbations outside the hori- 
zon until the universe passes through a phase of single 
field inflation that drives the isocurvature to zero. We 
should also stress at this juncture that a phase of sin- 
gle field inflation is not the only way to force pertur- 
bations into an adiabatic mode. In models where the 
non-adiabatic fluctuations persist through the end of in- 
flation, the fluctuations may become adiabatic during 
a phase of local thermal equilibrium with no non-zero 
conserved quantum numbers psl l46j . There are also 
other ways to achieve adiabaticity, such as the curva- 
ton scenario ^13, 16]. It is reasonable to suggest that in 
these models non-gaussianities generated from superhori- 
zon evolution of the curvature perturbation may be ob- 
servable. While we hope to explore these possibilities in 
future work, we believe that without understanding the 
complete evolution of perturbations until they become 
adiabatic, it is somewhat premature to claim observable 
non-gaussianity as a prediction of multifield inflation. 

The structure of the paper is as follows: We begin 
in section [ll] with a discussion of the SN formalism and 
the details of the model that we consider. In section 
mil we give the results for the two-, three-, and four- 
point statistics. In section IIVI we impose the condition 
that non-adiabatic fluctuations are damped away by a 
phase of effectively single field inflation, and we follow 
the evolution of tnl and ^nl during this process. We 
discuss the generalization of our work to n-point func- 
tions in section |Vl We conclude in section IVll Appendix 
El gives the details for the calculation of three and four- 
point statistics, and appendix IB] gives explicit expressions 
for non-linearity parameters describing flve and six-point 
functions. 



II. THE MODEL 

We use the same model as in [Ij] , and, as in that work, 
we closely follow the treatment of Vernizzi and Wands 
(irj (see also in our use of the SN formalism. What 
follows is a somewhat condensed version of the deriva- 
tions in [l| - the interested reader should see that work 
for details. 



We work with canonical two field models of the form: 



S = d X\/—g 
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a) 

Using (j) and x for the two scalar fields, the slow roll 
equations of motion are (assuming |ij| ^ and If/ial <C 
H\6J): 



3m2 



W . 



(2) 



The validity of slow roll throughout superhorizon evolu- 
tion is open to question. However, as argued in [ij, brief 
violations of slow roll should not substantially change the 
results. 

Because of their relative tractability, models of this 
type and their associated non-gaussianities have been 
well studied in the literature, for specific classes of poten- 
tials - for example sum-separable potentials (e.g. f47|). 
product-separable potentials (e.g. [H, |4^) and many 
others (e.g. [13, EJUl HI ) • 

As alluded to above, studying such a system is most 
readily done using the SN formalism jSOnSS] to calcu- 
late the evolution of the curvatu re p erturbation on uni- 
form density hypersurfaces, C, |54|. which is conserved 
for adiabatic modes whose physical size is greater than 
that of the horizon. The SN formalism relates C, at some 
comoving time tc to the perturbation in the number of 
e- foldings from an initially flat hypersurface aX t = t<t, 
generally taken to be horizon exit: 

C{t^,x) SN {t^,t,,x) = M {tc,t,,x) - N {tc,U) . (3) 

N is the unperturbed number of e-foldings, given by in- 
tegrating H from to tc'- 



N 



Hdt. 



(4) 



Viewing the number of e-foldings as a function of the field 
configuration on the hypersurface defined by horizon exit, 
0^(i*,x), and of tc, the perturbation in N can then be 
expressed in terms of the fluctuations of the scalar flelds 
at horizon exit, shown here to third order: 



SN : 



IJ 



NijScjj'Mi 



UK 



N 



ukc 



(5) 



The derivatives here are taken with respect to the flelds 
att — {Nj = -^jr) and SN is independent of the initial 
velocities since the potential is assumed to be slow roll 
at horizon exit. 



Statistics 

From the above formalism, one can straightforwardly 
derive expressions for cosmological observables in terms 
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of derivatives of the number of e-foldings. As in [l[ we fol- 
low the treatment of Vernizzi and Wands [13] and in ad- 
dition we also use the work of Byrnes, Sasaki and Wands 
(ioj - for a more detailed exposition of the derivation 
of SN statistics, please see those works. We begin by 
defining the curvature and scalar power spectra, V( and 



The spectral index is given by (the approximate equality 
denotes lowest order in slow roll): 



1 dlnT^c 



2 Eij'i'jNjjNj 



(8) 



e = —H/H^ is the usual slow roll parameter. Using the 
slow roll equations of motion, this can be written as 'ss'] : 



(6) 



nt -1 = -2e 



2ml WjjNjN, 



J 



2 
K 



(9) 

Three point statistics can be obtained in much the 
same fashion. The curvature bispectrum, _B^, is defined 
through: 



(10) 



Then, from ([3]), © and (HJ, we have: 



The bispectrum can be used to define the non-linearity 
parameter /nl [Hi which is the quantity most often ref- 
erenced in observational constraints: 



6. 



c 



(11) 



(7) From © and ©: 



UK 



IJKL 

r 



The star denotes a convolution and higher order terms 
have been neglected [^ \5% 

With some work (see [l[ for a more detailed exposi- 
tion) : 

i?(fci,fc2,fc3) = 47r47'2^i|| 



1 



\ 

with: 



^ , T.i.jN,iN^jNjj 



- ^ kik 



From the above and ((TT|) . /nl is: 



/NL 



(1 + /) + 



(13) 



(14) 



(15) 



The function / —1 — y^-p-j takes values between 

and 5/6 and depends upon the shape of the bispectrum 
[ll| . The first term of dTSl) is proportional to the tensor- 
to-scalar ratio, r, and so is guaranteed to be small [2l, 
thus only the second term can give rise to a large /nl^ 



5/NL- (^^iV2)2 



(16) 



This term, fj^^, is momentum-independent and local in 
real space, so it contributes to /^l^'j which, as mentioned 
in the introduction, can be a sharp probe of the number of 
dynamical degrees of freedom during inflation. However, 
as we have shown in [D], generating and keeping a large 
/nl^' along with an adiabatic power spectrum is far from 
easy. 

The leading contribution to the 4-point function is 
given by [ioj : 



4 



+ -N,A,A,N^B,B2N,cN^D 



kA2 



(3 perms) 



kJ't>kMJK) + i^ p™) 



(17) 



The subscript c denotes the connected part of the trispec- with 
trum; there is also a disconnected contribution (consist- 
ing of parahelograms of wavevectors) present even for 
gaussian fields. The trispectrum, T<^(fei, A;2, fes, fc4) can 
then be defined as: 

Although the general expression for is quite in- 
volved, it simplifies considerably in the case where we 
limit ourselves to models with independent gaussian fiuc- 
tuations at horizon exit (40| : 

rc(fei,fc2,fe3,fc4) = TNL [Pc{h3)Pc{h)Pdk4) + (11 perms)] 
54 

r5NL [Pc{k2)Pdk3)Pdk4) + (3 perms)] 

(19) 



ai = 



bi = 



EjNjjN. 



Nr 



3/2 ' 



1/2 



(23) 



Current experimental bounds are consistent with this in- 
equality, but should it be violated in future measure- 
ments, a large class of multifield models would be ruled 
out. Of course, such an inequality would be of less utility 
if the approach to adiabaticity wipes out any possibility 
of an observable (local) non-gaussian signal. 



\ki + kj\, and tnl and ^nl are given by: 



TNL 



5NL 



25 EijK NjjkNiNjNk 



54 



(20) 



The Potential 

Requiring that the SN formalism be of practical use 
places restrictions on the form of the potential, discussed 
in detail in [1] . The principal point is that in order to take 
derivatives of N with respect to the initial field values we 
need a way to easily relate the final and initial field values. 
This can be done if a suitable constant of motion can be 
constructed, which restricts the form of the potential to: 



W{<l>,x) = F{U{cl>)+V{x)), 



(24) 



Note that the parametrization is such that the coefficients 
multiply functions with different fc-dependence, making 
them observationally distinguishable. Restricting oneself 
to the above form gives the local-form non-gaussianity of 
the trispectrum. 

Before moving on we should point out that much of the 
recent interest in the trispectrum is due to the discovery 
of a useful inequality by Suyama and Yamaguchi, relating 
f^t to TNL [11 [ii: 



and also requires that there is a one-to-one mapping be- 
tween (j) and X- With W of this form, the appropriate 
constant of motion is: 



C 



1 



U' 



1 



y'ix 



ndx'- (25) 



The slow roll equations of motion are: 



TNL > 




(21) 



SHx'. 



F'U', 
F'V , 
1 



This is derived from the Cauchy-Schwarz inequality 



3TOp 



W{^,X) 



(26) 



(22) 



Here, as below, a prime denotes the derivative with 
respect the argument of the function. The slow 



roll parameters — {m'^/2){dW / d(l)iY and rj 



ij 
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mp{d'^W/d(t)id(j)j) are given by (we drop the repeated 
index on the diagonal 77 for conciseness): 



F 



(27) 



F J ' F'^ 
With the potential given above, N is now of the form: 

1 rwi^,xm^^ 



N 



2ml J, M^0(0,X(0))' 



1 



2ml J, WxWx),x) 



dx- (28) 



We impose a further restriction by requiring that the 
integral defining N can be divided into separate integrals 
over (j) and X: without using the functions (t){x) or x{4')- 
This gives two possibilities for the potential homogeneous 

W^[Ui^)+V{x)]\ (29) 

and exponential 

= WoExp[[/(0) + -t^(x)] . (30) 
For the homogeneous case, N is then given by: 



And for the exponential case: 

1 r 1 1 r 1 

N^~^ I ^^^d<^-^ / rrr— dx. (32) 



2ml J, U'{ct>) 



2mi 



V (x) 



When we end inflation we will also be interested in the 
parameters rj'^'^ and t^***, where a and s refer the direc- 
tions parallel and perpendicular to the inflaton motion in 
field space. These parameters control the masses of the 
adiabatic and isocurvature perturbation, respectively. In 
section llVl we will be interested in the case where the uni- 
verse passes through a phase of single field inflation where 
the isocurvature perturbations become heavy while the 
adiabatic ones stay light. These parameters are given by 



e 

_ e^-q't' - 2Vf*tXj^'t>x _|_ g0^x 



(33) 



Before moving on to our results and discussion thereof, 
we note that although our calculations assume slow roll, 
small deviations are possible. In particular, violations 
of slow roll are consistent with our 5N analysis so long 
as either the violation is sufficiently short or if is 
only weakly-dependent on the initial field value (attrac- 
tive trajectory) and friction is subdominant during any 
non-slow roll regime. The details of this argument can 
be found in There is also some numerical evidence 
to suggest the applicability of 5N beyond the slow roll 
regime - see [i^l and [Glj . 



III. RESULTS 

We include a summary of the results for both the ho- 
mogeneous and exponential potentials here. The details 
of the calculation are presented in appendix A of [H . We 
reproduce that exposition and expand it to include four- 
point statistics in appendix \K\ below. 



Homogeneous Potential: W{(j),x) ~ [U{4>) + V(x)T' 

Results here are expressed in terms of derivatives of the 
potential and the slow roll parameters defined in (l?fl) . As 
before a prime denotes a derivative with respect to the 
argument of the function - (j) for U and x for ^ ■ First, 
the mass parameters for the adiabatic and isocurvature 
fluctuations are given by: 



V 



(34) 



For the other quantities of interest, we first define: 
1 



Xh 



Vh 



1 



Vc4 - U,ei 



Vc4 - U,e} 



J7* -f K \ er. 
Then the observable T'(^ is given by 



and — 1 by 

— 1 = — 2e* — 



W* (xl , yl 



4 (xl , vl 



(35) 



(36) 



7 Ve' 



Vh 



2ei 



(37) 
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The full expressions for /nl, '''nl and ^nl are given 
in appendix]^ Many of the terms in these expressions, 
however, are multiplied by slow roll parameters. If we 
keep only those parts that are leading order in slow roll, 
we have (here and below we use e to refer to a generic 



first order slow roll parameter): 
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/ {Ua+Vgf (xH. _ ynV £i£i (Till 



(38) 




50l(t-t) (f/e + K)^^? 
5nl ^ {£*) + 7;z-- 



3 A 



7 



7(7-1) 
( 

7 



g-r jyA, 



7(7-1) 
( 

7 



2(7-1) ,. _ . ^ 2(7-1) 



(40) 



The ^ are defined by: 



I 4 7^ ([/ + y )' U'U'" + 372 (7 - 1) (iJ + 1/) C/'2[/" + 72 (7 - 1) (7 - 2) U" 



X2^ 4W^xExM _ 4 



4* = TO 



M/2 



7^ (?7 + Vf VV" + 37^ (7 -1){U + V) V'^V" + 7' (7 - 1) (7 - 2) 

([/ + y)^ 



(41) 



Exponential Potential: W{(j), x) = WoExp [(7(0) + V"(x)] 



define: 



The slow roll parameters (j27p take a particularly sim- 
ple form for an exponential potential of the form pop . 
and readily lead to the following expressions for the mass 
parameters yy*'* and 77°''^: 



2et 2e} 

— , ye = — 

Then the observables Vq and 77,^ — 1 are given by: 



— • (43) 



(44) 



e'^7?'^ 



4e0gX + ^XrjX 



£Xj^0 - 4e*eX + e'^??^ 



(42) 



71^ - 1 = -2e, - 4 



(45) 



Next, in a similar fashion to the previous section, we 



As before full expressions for /nl, tnl and ^nl can be 
found in appendix but if we keep only the leading 
terms in slow roll, we have: 



7 



6 .(4) 



yi 

2tl 



(46) 



64 



TNL 0(£*) + 



(47) 



.gNL O (e 



200 l^e^ 



( - 



(48 ^ - 26ejr7,x + 2 ry,^ ^ ^ Q ') + ("48 ^ + Gry,^,^ - 4 ry,^ ^ ^ 22e? {^4 + + ') 



(48) 



IV. DAMPING AWAY ISOCURVATURE AND 
THE FATE OF NON-GAUSSIANITIES 



Homogeneous Potential: W{(j),x) ~ [U{4>) + V'(x)]'' 

We are now in a position to assess the magnitude of 
(7nl and tnl as isocurvature is damped away, beginning 
with the homogeneous potential. We are only interested 
in the evolution of non-slow roll suppressed parts of the 
local trispectrum, so we will focus on the leading terms 
given in §^ and 

In order to translate the correlation functions of pri- 
mordial fluctuations into observations, we must follow 
their evolution until they become conserved quantities, or 
until they are observed. The latter would require detailed 
knowledge of the cosmological history from the present 
all the way back to the inflationary epoch. Given our 
ignorance of the details of much of the early universe, we 
focus on the former approach. Accordingly we wish to 
study TNL and ^nl as they are forced into constant val- 
ues and, since correlations of the curvature perturbation 
will evolve outside the horizon as long as non-adiabatic 
fluctuations are present, this means we need to evalu- 
ate (|39p and (|40|) after any non-adiabaticity has been 
damped away and only the adiabatic mode is left. Here 
we achieve this by requiring that there is a short phase of 
effectively single field inflation sometime before reheating 
while the observationally relevant modes are outside the 
horizon. We repeat that this is not the only possibility 
for achieving adiabaticity, but this choice ensures that 
we are able to calculate observables without making any 
assumptions about post-inflationary dynamics. In [l| we 

perform the same analysis and find that becomes 
slow roll suppressed as we approach adiabaticity. 

We now consider the equation of motion for the isocur- 
vature fluctuations which at leading order in slow roll 



takes the form [GC 



W 

5s + 3H6s + -^Tj^'Ss + —6s ^ . 



The solution of this equation is [60|, |62 



<5scxa(t)-3/2 



aH 



where v is given by 



1.2 = - 



9 Wri" 



4 jripH'^ 



(49) 



(50) 



(51) 



So for 77"^ > |, we find 



|<5s| cx a(t)-3/2 
\Ss\ - Exp 



(52) 



and the isocurvature fluctuations are rapidly damped 
away. 

We will now examine the conditions for large 77'*'', de- 
fined in dMl). RecaU that ^ gives 



.4>x ^ 2— '^—\f^ 

7 



(53) 



and (from (1571) ') I7I must be 0(1) or larger to guarantee 
scale invariance. Thus when 



(IMD 



is large, we must have either f]'^ }^ or tj^ > . If both 
r]'^ and rj^ are large, then rj"'"' will also be large, and 
inflation will quickly end before the non-adiabatic modes 
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have been damped away, so we will not be interested in 
this case. 

Now let us consider tnl, which we write in the follow- 
ing simplified form: 



(54) 



Given the appearance of 77*'', one might expect that for 
sufficiently large 77'*'* (and large yy''^ is required to damp 
away isocurvature) observably large tnl could be pro- 
duced even as isocurvature is damped away. However, 
the coefficient of 77^***, specifically is a dynamical 

quantity whose time-dependence is affected by the be- 
havior of T]"'^ . To find this time-dependence, we first note 
that: 

— = -2e'*77^ - 2V^7y*x ^_ 4(g0)2 _^ (55) 
H 



— = -2e^T]^ - 2V^ri'^^ + 4:{e^)^ + 4e'^eX. (56) 
H 

Now, for the damping of isocurvature either r;"^ > 1 or 
ij^ > 1. In the former case we can neglect all but the 
first term on the right hand side of ([55]) . and so we find 
that 



e'^(i) cx Exp 



(57) 



Similar remarks apply to e^, in the case that t]^ > 1. 

It it is then straightforward to write down the evolution 
of Tnl as isocurvature is damped away: 



rNL-0(e*) + e^(l) X (77'') Exp 



4 / C^Hri'^dt 



(58) 



Cri is a number which is always greater than 1 whose 
value depends on the particular direction of the effective 
single field during this phase. We thus conclude that 
Tnl will always be slow roll suppressed upon entering 
the purely adiabatic solution after a phase of effectively 
single field inflation. 

The argument is similar for ^nl- By examining (1401) . 
we can see that the part of (/nl which is not slow roll 
suppressed is proportional to ^^f-- As we argued above, 
either e*^ or will be exponentially damped when 77"* 
becomes large. In order for (;nl to remain large as non- 
adiabatic fluctuations are damped away, one of the terms 
in square brackets in (j40l) must grow exponentially. Since 
we must have e < 1 throughout inflation, we cannot have 
exponential growth from terms which contain only e*^ and 
e^. An exponentially growing r]'^ would result in a double 
exponential damping of e"^, as can be seen from (1571) . and 
thus a rapid damping of ^nl (and likewise for 77^). 

This leaves only the terms proportional to How- 
ever, the evolution of is not independent of that of the 



other slow roll parameters. Specifically, we note that 



n 



+2e'''r]^ + -e^rj'f' + -V^r]'!'^ (59) 
H 7 7 



H 



7 



7 



(60) 



If i^"^ were positive and exponentially growing, this 
would mean that rj'^ would quickly become negative and 
exponentially growing. This would lead to a double ex- 
ponential growth of e"^ which would result in a quick 
end to inflation before the non-adiabatic fluctuations are 
damped away. If were negative and exponentially 
growing, then rj'^ would exponentially increase as well, 
which in turn would result in a double exponential sup- 
pression of e"^, and thus a rapid damping of (?nl- The 
same analysis holds for So we find that like f'^^ 

and Tnl, .9nl will always be slow roll suppressed upon 
entering a purely adiabatic solution after a phase of ef- 
fectively single field inflation. 

We pause to note that while the terms labeled ©(e*) 
are proportional to slow roll parameters, there are cases 
when they are not negligible. As shown in Appendix 
[A1 these terms contain several factors which depend on 
the details of the potential. In fact, there are examples 
showing that these terms can be significant for f'^^ (63l . 
[63 |. and a similar mechanism should apply for higher 
point statistics as well. 



Exponential Potential: W{(j),x) = WoExp[U{(l>) + V{x)] 

An entirely analogous damping takes place in the case 
of the exponential potential. Once again the dominant 
contribution to tnl will be proportional to yy^'' ^ 

rNL^O(e*) + 0(l)x r,^' , (61) 



and 77*'' is 



A large 77** will then imply (again) either a large 77"^ or 
a large 77'"^ and lead (in the same fashion as for the ho- 
mogeneous potential) to the exponential damping of the 
leading term in tnl as in ([58]) . 

The argument for (^nl is similar. The exponentially 
suppressed factor also multiplies the leading order 
term in (Jnl and the factors in the expression cannot be 
exponentially large (by the same arguments as above). 
Thus the leading term in ^nl is driven to zero as the 
fluctuations are forced into an adiabatic mode. 



V. HIGHER POINT FUNCTIONS 

Although we have examined only the three- and four- 
point functions in detail, we can extend our analysis to 
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higher point functions as welL The n-point function is 
always proportional to n factors of N with at least one 
derivative acting on each N. If we restrict our attention 
to models where the scalar field fluctuations are gaus- 
sian and independent at horizon exit, then the contribu- 
tions to the non-gaussian correlation functions are pro- 
portional to the various ways of contracting the deriva- 
tives of N with Kronecker deltas. If we focus on the con- 
nected part of the non-gaussian correlations, then there 
is no subset of the factors of N whose derivatives are 
contracted independently of the other factors. In this 
case the leading contributions to the n-point function 
will contain 2{n — 1) derivatives. We will briefly discuss 
the sub-leading contributions below. 

For a general n-point function we can characterize the 
various local forms by defining a set of non-linearity pa- 
rameters. The number of independent non-linearity pa- 
rameters is given by the number of distinct ways of apply- 
ing 2(n — 1) derivatives to n factors of N with at least one 
derivative acting on each factor of N. This is equivalent 
to the number of "free trees" constructed from contract- 
ing n vertices with n — 1 edges, called t„ in graph theory 
and combinatorial analysis |65l . |66| . We can write the set 
of non-linearity parameters for the n-point function as 



p(n) _ 



E 



N 



A1A2. 



N 



A1A3 



Ar2 



(62) 

where the index i runs from 1 to t„ and labels the various 
ways of distributing 2(n— 1) derivatives in the numerator. 

We can make the identifications i^A 

54 , 



_ 6 f (4) p(4) 



TNL, and i^L 2 



NL,1 — 5-' 

2^ffNL- There are 3 parameters for the 
five-point function, 6 for the six-point function, 11 for 
the seven-point function, 23 for the eight-point function, 
and so on [gl, [6^ . We give explicit expressions for the 
five and six-point non-linearity parameters in appendix 

IE 

The structure of the derivatives of N follows a pat- 
tern which allows for an easy determination of the most 
important contributions to each of the non-linearity pa- 
rameters. Specifically, the mth derivative of N will be of 
the form 



term which is lowest order in slow roll parameters will 



For a product of n factors of N 

with a total of 2(n — 1) derivatives, the leading term is 
O (^e7^"~^^^ , and so the leading term is the only one 
which is not automatically slow roll suppressed. 

The leading term for all second and higher derivatives 
of iV contains a factor -£-5^ . One can see this directly for 
the second and third derivatives by examining equations 
(|A10IIA14|) for the homogeneous potential and equations 
(jA20IIA24|) for the exponential potential. It is straight- 

forward to verify that derivatives of with respect to 
0* and X* are proportional to a sum of slow roll param- 
eters multiplied by the combination It is exactly 

this quantity which becomes exponentially suppressed as 
non-adiabatic fluctuations are damped away by passing 
through a phase of single held inflation. As long as the 

terms which are multiplying are not exponentially 
growing during the phase of single fleld inflation, we flnd 
that the leading term in all second and higher derivatives 
of N becomes exponentially suppressed as non-adiabatic 
fluctuations are damped away. As we discussed above, 
the various slow roll parameters do not evolve indepen- 
dently, and thus generally cannot grow exponentially 
without affecting the evolution of e"^ and e^. As a re- 
sult, we can conclude that the leading contribution to all 
local form non-gaussian n-point functions becomes sup- 
pressed as non-adiabatic fluctuations are damped away 

fn) 

during a phase of single field inflation, and thus all F^-^ ^ 
will be slow roll suppressed. 

There is a complication which we have so far over- 
looked, which is that there are other contributions to the 
non-linearity parameters which contain more than the 
minimum number of derivatives acting on the factors of 
N, see for example 63, 67 69j. These terms are generally 
known as loop corrections because they involve integrals 
over internal momenta. These terms all contain second 
and higher derivatives of N, and are thus suppressed for 
the same reason as the leading terms discussed above. 



fc=0 



[si j 



(63) 



Each of the non-linearity parameters describing the n- 
point function has n factors of N with a total of 2 (n — 1) 
derivatives in the numerator and 2(n — 1) factors of 
N^K in the denominator. This means that the denom- 
inator will be O ^e, so only the terms which 

are O ^e* ^" or larger in the numerator will make 

a contribution which is not automatically slow roll sup- 
pressed. Now if we have a product of n factors of N 
of the form M"^^^ N^""^^ ■ ■ ■ N^"'"^ where N^""^ refers 
to the TTiith derivative of N with respect to (j>l, the 



VI. CONCLUSIONS 

We have extended the analysis of [l| to include the 
calculation of the trispectrum. We calculated tnl and 
(7nl for inflationary potentials of the form W{4>, x) = 
[C^(0) + Vix)r and W{q^,x) = Wo^xp [£/(</.) + Vix)] by 
using the SN formalism. We focused on the case in which 
there is a phase of effectively single fleld inflation which 
damps away non-adiabatic fluctuations before the end 
of inflation. This ensures that all of the n-point statis- 
tics are conserved during the subsequent evolution of the 
universe, and it also guarantees that the power spectrum 
is purely adiabatic as indicated by observation. Under 
these conditions, we flnd that tnl and (/nl are always 
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slow roll suppressed upon entering a purely adiabatic so- 
lution. We also discussed the extension to higher point 
functions and argued that all of the non-linearity param- 

(n) 

eters i^L i ^^^^ describe the local form non-gaussianity 
for any n-point function will be slow roll suppressed after 
a phase of effectively single field inflation. 

There are, of course, some serious limitations with our 
result. In particular, one might worry about different 
potentials, different (slow roll violating) trajectories and 
different approaches to adiabaticity. As discussed in [l|, 
there are some arguments which suggest that the first 
two of these would not have significant effects on our 
results. In yy] (and as recapped at the end of section HI)) 
we show that small violations of slow roll do not affect the 
results and argue that larger ones would be inconsistent 
with the observed power spectrum. With regards to a 
different potential, it is unclear how such a change could 
alter the fundamentals of our analysis, though without a 
detailed study one cannot be certain of this. 

It is, however, the particular approach to adiabatic- 
ity that is the most restrictive assumption we make. 
While our results illustrate the challenges associated with 
generating an adiabatic spectrum and local form non- 
gaussianities, there a number of scenarios where the fluc- 
tuations in the cosmological fluid enter an adiabatic mode 
without a phase of single field inflation. As discussed in 
the introduction these include (amongst others) curvaton 
models and a phase of thermal equilibrium. We hope to 
explore the evolution of the bispectrum, trispectrum and 
beyond in such scenarios in forthcoming work. 
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Appendix A: Details for Homogeneous and 
Exponential Potentials 



Varying N then gives: 
1 



dA^ = 



1 



u'j ^ d(f>, \u'j^ d<f>, \v' 



(Al) 



Note that in deriving the above we had to account for 
the dependence of 0c and Xc on both 0* and x* ■ We will 
also need: 



d0e / dC . 

_dxcfdC_ 

dc [deb.' 



dC , 

dC , 
a — dx* 

ox* 



(A2) 



From (|25|) we have: 

dC _ 1 



p 
UL 



dC ml 



1^ = ^^ (A3) 



and 



dC 



U' 



dC _ ml 



The time t(. defines a surface of constant energy: 

W {4'c, Xc) = constant . 
Differentiating with respect to C then gives: 



dC 



(A4) 



(A5) 



(A6) 



Using the above and Wtt,/Wx = U' /V , we can differen- 
tiate the expression for C in and obtain (after some 
manipulation) the following expressions: 



dC 

dXc 



mt 



dC ml 



c/'( — + — 

, / 1 1 

V \ 



1 uiv. 



ly'2 

c* c 



ml c/;2 + 

ml + ■ 



(A7) 

Substituting (|A3|) and (|A7j) into ()A2p allows us to read 
off the following: 



Homogeneous Potential: VK((^,x) = [f/(<^) + V{(l>)]'^ 
With N and C given by (^5]) and ([5T |) : 



N 



1 



1 



"/ml J. U' (0) 7to2 V (x) 



'2 



__ _ U', d(j)c 

W*^ WTWul' dx* 

dxc U7 Vc dxc_ 

8(1), u;? + K'' ' dx* 

The derivatives of N are then: 

dN _ 1 Xh dN 
dx* 



v. 



'2 



UL 



c/;2 + v!^ v! ■ 



1 Vh 



(A8) 



(A9) 
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We have used the slow roll parameters defined in (1271) : 
while Xfi and yh are defined in ([35]): 



1 



Xh 



Vh 



1 



In a similar vein, we can find the second derivatives: 
927V 1 



7m2 



1 - Xh 



2et 



7 + 1 



S^iV 1 



7m2 



1 f '^'^* ,1 



7m2 



- 1 
(AlO) 



and third derivatives: 



d^N _ 1 
7m2 



1 r;? ,7-1 



p \/2ef 7*™? 



Xh / 7 



7 V* ,7~1 ^*_,7~1, 



26*" 7"ip 



1 77,^ 
'"p ^2^ 



TOp \^ 2 V2 g0 72g0 3/2 



1 (C/e + Kf efe? [2 / .2 



A , 4(7-1) 



"^f V2ef^'(C^* + ^*f 1 7 



7 



(All) 



7m2 



1 77? ^ 
^2^ jmp 



'2d 



Vh 



7 



X 2 
7 ^* 



l2A/2eX3/2 72eX3/2 



7 - 1 



^2ei 



7-1 
irup 



'2d 



\3 



'2ei 



+ 



7 



"^P V2ej'/' (t/* + K)' 
+ 2 %x 2 _ 2 ,x 2 ^0 2 ^ 14(7 -l) ^xg0 (,x^0 _ ,0^x) + ,0 2 2 _ 2 ^0 

6£Xg*^«s ( 77"^ - 77X + 



gX _ 



(A12) 
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1^1 



- 1 



7 1 (C/c + V<,f efe^ ( 2 



-el e-^ - 2 r;-^ I + 



4( 7-1) 
7^ 



.2(7-1) 



_ 



(A13) 



dxld(t)* iml 



1 ?7? 7 - 1 



'2d 



{U, + y,)' e^A 1 



7 1 (C/c + V,f e^et r 2 



7^70 



3 .4 



gX 2 ^0 _ g0 ^x 



4(7-1) 
7^ 



(A14) 



where ^ is defined by (|1T|) : 



^2^ 4W^^W^_ 47' (^ + + 37' (7 - 1) (£/ + t^) C/^'C/" + 7' (7 - 1) (7 - 2) U'^ 



= m'. 



and 77^" is given in (|34 



1^2 



TO, 



((7 + y)' 



4 " ,x ^"^ ,xxx 



T4/2 



4 7^ {u + y)" rr" + 37' (7 - 1) ([/ + T^) V'^V" + 1) (7 - 2) y 

(c/ + y)' 



V 



(El 



From these it is straightforward to use the relevant SN equations to obtain expressions for /nl, tnl and g^h- 
( 



6 (4)_ 2 

5^NL - ^ 



1 



( 



- 7 + 1 .T/i 



Vh 



1 - 



2ei 



- 7 + 1 y/i 



V 



/ ^1 I vl 



( VhS^ tee? ( inT 

(l + S)' 



(A15) 
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and 



.9nl 



50 1 
277 



yl 



+ 3 



24 



7-1 

7 

^7^ 

7 



7^2 



7 



(7 - 1) Tyf , 7-1 



A el 



Vh 



7 ?7: 



2d 

X2 



2e^ 



7 



(7-1)77? , 7-1 



4e? 



2 e? 



2ei 



( Vh^h _ fuL^ 7-1 

[etei eiVW 7 



4''' „x 



^ll_ilg0gx (gx _ g0) + gx 2 ^0 

(gX2^ 



g0 jjX ^j^X 



4>^ _ ^4>^ ^x-2 



7(7-1) 
f 

7 

2(7-1),. , 2(7-1) 



7(7- 



(A17) 



Eliminating the slow roll suppressed terms gives the expressions in ([55)1 . ([5^ and (HUl) . 



Exponential Potential: VK((^,x) = WoExp [(7(0) + V{x)] 
With iV from ^ and C again given by ([25]), 



00 



N 



1 



1 



1 



1 



2ml J, ^'W 2m2 7, ^ (x) 

the variation gives: 
1 



dx, (ESI) 



dN 



2m2 



1 



2m2 



(A18) 



The analysis begins in the same fashion as in the previous 
section, with the expressions for the derivatives of (pc and 
Xc with respect to 0, and x*j given in (jA8l) . From this 
and (|A18I) we obtain the following expressions for the 
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first derivatives of N: 

dN 1 Xe 



2m 



P \/2€l 



dN 
dx* 



1 



Ve 



(A19) 



The slow roll parameters are defined in ([27]). Xe and ye 
are defined in (1431) as: 



Xe 



24 _ 2£X 

The second derivatives are then given by: 



d'^N 
dxl 

d(t>*x* 



1 1 
2™^ 2ei 
1 1 



2ml 2e? 



2ef Xe 



(??? - 2eJ)2/e 



4g'^gX 

,2 Ic 



4g<l'gX 
.2 'c 



2ml ^, 



1 2efe} 

'Ic 



6j 6^ ^1 



(A20) 



And the third derivatives are: 
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2e, 



(2 (et - rjt) vt + ) + 3 (2.t - r^f) ^ry; 



+ ef (48 ^ + 226? (r;,^ - ,7?) ~ 677,^? + 4 ,7? ^ _ 2^ _ ^0 3 2 _ ^g^^^^x + 2 ,7? ^ _^ ^ ') 



(A21) 



§- (2 (e? - Vi) + e + 3 (2eJ - r?? ) 



c c _ss 



(2 rj^ ' + e + (-26ef r?? - 4 77? ^ + 677?/;,^ + ^ 



+ 2 et (48 e,^ ' + 22ef (r;? - ,7,*) - 677^,7,^ + 4 77^ - ') - ^ (48 ' - 26e,^,^ + 2 77^ + 4^ ') 



(A22) 



f0pX 



(2 ?7?' + e 



^ ,0 2 (_26e^77? - 4 77? V 677^77,^ + ' 



+ ^ et (48 ' + 22ef (77? - 77,*) - 677,^77,^ + 4 7?,* ' - 4^ - (a8 ^ ' - + 2 77^ + 4^ 



, (A23) 
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'(26j~77j)^7;r + ^ 

6-^ 6-^ \ 



(2 4 ' + 4^ ') ^4 ' (-266^77,^ ~ 4 7;f + 677^? + S,t ') 



+ ef (48 ^ + 22e? (^f _ ^x) _ 677^77? + 4 77,^ ^ - Q ^) - ef ' (48 ^ - 26e^vi + + ') 

77'"' is given in (|42]) : 



£X^0 - 4e'^e>^ + e'^77>^ 



As above, the SN equations then give the observables of interest 

2 



6 Ai) _ 



3 ( ')f-2e: » 



{24 ^ 2^ 



(A24) 



(A25) 
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+ 8^7;- {y^et - Xeej) (x^ ^ {r,t - 2et) + yl ef{r,i - 2e^)) 



(A26) 
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and 



100 



5nl = 



[yiei + xid 



(2 (et - vi) 4 + it ') + Vl^l 4 ' (2 (ej - ili) + e ')) 



+ et ' (48 ^ - 26eM + 2 + ^ ') + ("48 ' + Gvtvi " 4 J??' + 22e? (-r;,^ + ry?) + 



(A27) 



Suppressing the slow- roll terms gives the expressions in p6)) . (|47l) and (|48 
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Appendix B: Five- and Six-Point Non-Linearity 
Parameters 



In this appendix we give explicit expressions for the 
non-linearity parameters which describe the amplitude of 
the various local limits of the five- and six-point functions 
assuming independent gaussian scalar field fluctuations 
at horizon exit. By using the general expression (|62p we 
obtain for the five-point non-linearity parameters: 



And for the six-point non-linearity parameters: 



r^(6) 
NL,1 



p(6) 
NL,2 
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NL,3 
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M 
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ElJKLM^JJKLMNjNjNjjNxNM 



(B2) 



It can easily be shown that there are many other con- 
sistency relations like the one found by Suyama and Yam- 
aguchi [58, 59] which relate the non-linearity parameters 
for higher point functions. For example, if we take 
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NL,1 
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then the Cauchy-Schwarz inequality 
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